Multidynamic Crystalline Molecular Rotors Comprising an N-Heterocyclic Carbene Binuclear Au(I) Complex Bearing Multiple Rotators

We report multidynamic molecular rotations in crystals using a concave-shape N-heterocyclic carbene (NHC) binuclear Au(I) complex rotor bearing pyrazine and tetrahydrofuran (THF) molecules as multicomponent rotators. Single-crystal X-ray diffraction (XRD) measurements revealed that two THF molecules are located near the central pyrazine encapsulated by two bulky NHC ligands. From ²H solid-state NMR analysis, it was observed that the pyrazine rotated in a 2-fold site exchange with a 180° rotational angle and a 31 kJ mol⁻¹ energy barrier, while the THF molecules showed a 23°-38° libration with a lower energy barrier (14 kJ mol⁻¹). Interestingly, the pyrazine rotation was accelerated when the THF molecules rotated in fast site exchange with a large angle of libration, suggesting that the rotators exhibit multidynamics in a correlated manner.     

Probing Electronic Symmetry Reduction during Charge Migration via Time-Resolved X-Ray Diffraction

The present work focuses on probing ultrafast charge migration after symmetry-breaking excitation using ultrashort laser pulses. LiCN is chosen as prototypical system because it can be oriented in the laboratory frame and it possesses optically-accessible charge transfer states at low energies. The charge migration is simulated within the hybrid time-dependent density functional theory/configuration interaction framework. Time-resolved electronic current densities and simulated time-resolved x-ray diffraction signals are used to unravel the mechanism of charge migration. Our simulations demonstrate that specific choices of laser polarization lead to a control over the symmetry of the induced charge migration. Moreover, time-resolved x-ray diffraction signals are shown to encode transient symmetry reduction at intermediate times.     

Ring-polymer Molecular Dynamics Simulation for the Adsorption of H2 on Ice Clusters (H2O)n (n=8, 10, and 12)

In the interstellar medium, the H₂ adsorption and desorption on the solid water ice are crucial for chemical and physical processes. We have recently investigated the probabilities of H₂ sticking on the (H₂O)₈ ice, which has quadrilateral surfaces. We have extended the previous work using classical MD and ring-polymer molecular dynamics (RPMD) simulations to the larger ice clusters, (H₂O)₁₀ and (H₂O)₁₂, which have pentagonal and hexagonal surfaces, respectively. The H₂ sticking probabilities decreased as the temperature increased for both cluster cases, whereas the cluster-size-independent profiles were observed. It is thought that the size independence of the probabilities is qualitatively understood from the similar binding energies for all the three cluster systems. Furthermore, the RPMD sticking probabilities are smaller than the classical ones because of the reduction in the binding energies owing to nuclear quantum effects, such as vibrational quantization.     

Understanding Conformational Properties and Role of Hydrogen Bonds in Glycosaminoglycans-Interleukin8 Complexes in Aqueous Medium by Molecular Dynamics Simulation

Atomistic molecular dynamics simulations were performed under ambient conditions to explore the conformational features and binding affinities of hexameric glycosaminoglycans (GAGs) with chemokine Interleukin8 (IL8) in an aqueous medium. We tried to understand the role of hydrogen bonds (HBs) involving conserved water in mediating the interactions. The Luzar-Chandler model was adopted to study the kinetics of HB breaking and formation concerning different water-mediated HBs. The conformational flexibilities of bound GAGs are due to the flexible glycosidic linkages than the occasional/rare ring pucker conformation. The free energy landscape constructed with ϕ, and ψ, depicted that different conformational minima associated with the glycosidic linkage flexibility of the GAGs in bound states are separated by energy barriers. The binding affinities of IL8 towards GAGs are favored through the electrostatic and non-polar solvation interactions. 4-different types of conserved water were explored in the solvent-mediated binding of GAGs with IL8. The average lifetime of the IL8-GAG direct HB pairs was ∼ten times less than the IL8-GAG-shared water HBs. This is due to the rapid establishment of HB breaking and reformation kinetics involving water of a shared layer. We find that despite the highly negatively charged surface of GAGs, the IL8 surface populated by non-cationic amino acids could serve as a promising binding site in addition to the cationic surface of the protein.     

Coupling Effects of Electrostatic Interactions and Salt Concentration Gradient in Polymer Translocation through a Nanopore: A Coarse-Grained Molecular Dynamics Simulations Study

We study the influence of polymer pore interactions and focus on the role played by the concentration gradient of salt in the translocation of polyelectrolytes (PE) through nanopores explicitly using coarse-grained Langevin dynamics simulations. The mean translocation time is calculated by varying the applied voltage, the pH, and the salt concentration gradient. Changing the pH can alter the electrostatic interaction between the protein pore and the polyelectrolyte chain. The polymer pore interaction is weakened by the increase in the strength of the externally applied electric field that drives translocation. Additionally, the screening effect of the salt can reduce the strong charge-charge repulsion between the PE beads which can make translocation faster. The simulation results show there can be antagonistic or synergistic coupling between the salt concentration-induced screening effect and the drift force originating from the salt concentration gradient thereby affecting the translocation time. Our simulation results are explained qualitatively with free energy calculations.     

Role of the Residue Q1919 in Increasing Kinase Activity of G2019S LRRK2 Kinase: A Computational Study

Mutations in multi-domain leucine-rich repeat kinase 2 (LRRK2) have been an interest to researchers as these mutations are associated with Parkinson's disease. G2019S mutation in LRRK2 kinase domain leads to the formation of additional hydrogen bonds by S2019 which results in stabilization of the active state of the kinase, thereby increasing kinase activity. Two additional hydrogen bonds of S2019 are reported separately. Here, a mechanistic picture of the formation of additional hydrogen bonds of S2019 with Q1919 (also with E1920) is presented using ‘active’ Roco4 kinase as a homology model and its relationship with the stabilization of the ‘active’ G2019S LRRK2 kinase. A conformational flipping of residue Q1919 was found which helped to form stable hydrogen bond with S2019 and made ‘active’ state more stable in G2019S LRRK2. Two different states were found within the ‘active’ kinase with respect to the conformational change (flipping) in Q1919. Two doubly-mutated systems, G2019S/Q1919A and G2019S/E1920 K, were studied separately to check the effect of Q1919 and E1920. For both cases, the stable S2 state was not formed, leading to a decrease in kinase activity. These results indicate that both the additional hydrogen bonds of S2019 (with Q1919 and E1920) are necessary to stabilize the active G2019S LRRK2.     

Rapid separation and binding configuration prediction of the components in Danshen decoction to endothelin A receptor using affinity chromatography and molecular dynamics simulation

As a famous traditional Chinese formula, Danshen Decoction has the potential to relieve the pain of pulmonary arterial hypertension patients, however, the functional components remain unknown. Herein, we reported a method to screen the functional components in Danshen Decoction targeting endothelin receptor A, an accepted target for the treatment of the disease. The receptor was functionalized on the macroporous silica gel through an epidermal growth factor receptor fusion tag and its covalent inhibitor. Using the affinity gel as the stationary phase, the bioactive compound was identified as salvianolic acid B by mass spectrometry. The binding kinetic parameter (dissociation rate constants k_d) of salvianolic acid B with the receptor was determined via peak profiling. Using the specific ligands of the receptor as probes, the binding configuration prediction of salvianolic acid B with the receptor was performed by molecular dynamics simulation. Our results indicated that salvianolic acid B is a potential bioactive compound in Danshen Decoction targeting the receptor. This work showed that receptor chromatography in combination with molecular dynamics simulation is applicable to predicting the binding kinetics and configuration of a ligand to a receptor, providing crucial insight for the rational design of drugs that recognize functional proteins.     

4-硝基咪唑A-带激发态结构动力学

通过共振拉曼光谱实验和量子化学计算的方法研究了4-硝基咪唑(4NI)A-带激发态衰变动力学. 对4NI的振动光谱、紫外电子吸收光谱、荧光光谱和共振拉曼光谱进行了指认. 在全活化空间自洽场法(CASSCF)/6-31G(d)计算水平下获得了单重激发态S₁(n_Oπ^*)和S₂(ππ^*)和势能面交叉点S₁(n_Oπ^*)/S₂(ππ^*)的优化几何结构和能量, 分析了A-带共振拉曼光谱的强度模式特征, 获得了短时结构动力学, 并结合全活化空间自洽场法(CASSCF)理论计算结果确定了4NI 在S₂(ππ^*)态衰变通道主要是S_{2, FC}→S_{2, min}(ππ^*)→S₀辐射弛豫.     

分子动力学模拟研究[Bmim][PF6]+水+醇三元体系的微观结构和分子间相互作用

研究离子液体体系的微观结构和分子间相互作用具有重要意义. 本文对1-丁基-3-甲基咪唑六氟磷酸盐([Bmim][PF₆])+水+乙醇和[Bmim][PF₆]+水+异丙醇三元体系进行了分子模拟研究, 计算了径向分布函数和不同组成的水-醇混合溶剂与离子液体阴阳离子间的相互作用能, 并将其分解为库仑相互作用能和Lennard-Jones(LJ)势能. 在此基础上, 研究了溶液体系的微观结构、分子间相互作用和相行为. 结果表明, 水倾向于与离子液体阴离子和阳离子极性部分作用, 醇倾向于与阴离子和阳离子非极性部分作用; 库仑力主导阴离子-溶剂相互作用, 色散力主导阳离子-溶剂相互作用, 阴阳离子的缔合状态对色散力影响较小, 对库仑力的影响非常显著.     

Structure and Dynamics of Oligonucleotides in the Gas Phase

By combining ion‐mobility mass spectrometry experiments with sub‐millisecond classical and ab initio molecular dynamics we fully characterized, for the first time, the dynamic ensemble of a model nucleic acid in the gas phase under electrospray ionization conditions. The studied oligonucleotide unfolds upon vaporization, loses memory of the solution structure, and explores true gas‐phase conformational space. Contrary to our original expectations, the oligonucleotide shows very rich dynamics in three different timescales (multi‐picosecond, nanosecond, and sub‐millisecond). The shorter timescale dynamics has a quantum mechanical nature and leads to changes in the covalent structure, whereas the other two are of classical origin. Overall, this study suggests that a re‐evaluation on our view of the physics of nucleic acids upon vaporization is needed.